Spin-wave instabilities in spin-transfer-driven magnetization dynamics 
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We study the stability of magnetization precessions induced in spin-transfer devices by the injec- 
tion of spin-polarized electric currents. Instability conditions are derived by introducing a general- 
ized, far-from-equilibrium interpretation of spin-waves. It is shown that instabilities are generated 
by distinct groups of magnetostatically coupled spin-waves. Stability diagrams are constructed as a 
function of external magnetic field and injected spin-polarized current. These diagrams show that 
applying larger fields and currents has a stabilizing effect on magnetization precessions. Analytical 
results are compared with numerical simulations of spin-transfer-driven magnetization dynamics. 

PACS numbers: 75.60.Jk, 85.70Kh 



Currents of spin-polarized electrons can induce large- 
amplitude magnetization precessions at microwave fre- 
quencies in small-enough magnetic devices 

HQ. There 

is mounting experimental evidence that these so-called 
spin-transfer phenomena do occur in nano-pillar or nano- 
contact devices under current densities of the order of 
10 6 - 10 8 A/cm 2 This discovery has boosted the 

already widespread interest in the physics of the interplay 
between magnetism and electron transport, and has trig- 
gered efforts toward the promising development of new 
generations of microwave spin-transfer nano-oscillators. 

A spin-transfer device is a non-linear open system, 
driven far-from-equilibrium by the action of the spin- 
polarized electric current. The excited magnetization 
precessions represent strong excitations of the magnetic 
medium, which in principle may give rise to various types 
of instability and eventually to transitions to chaotic dy- 
namics. A parallel can be drawn with ferromagnetic- 
resonance Suhl's instabilities 7], in which certain spin- 
waves can get coupled to the uniform precession and start 
to grow to large non-thermal amplitudes, thus destroying 
the spatial uniformity of the original state. 

In this Letter, we demonstrate that spin-wave insta- 
bilities may occur in spin-transfer-driven magnetization 
dynamics as well. However, the system is far from equi- 
librium and the classical notion of spin-waves fails. In- 
deed, it is the large- amplitude magnetization precession 
induced by spin transfer that plays the role of reference 
state, and spin-waves only exist in a generalized, non- 
equilibrium sense, as small-amplitude perturbations of 
that state [il-fioj]. This scenario emerges with clarity in 
the time-dependent vector basis in which the reference 
magnetization precession is stationary. The spin-wave 
equations in this basis are characterized by two features: 
(i) a well-defined dispersion relation uj (q;cos8 ), whose 
non-equilibrium nature is revealed by its explicit depen- 
dence on the magnetization precession amplitude cos#o; 
and (ii) the presence of time-periodic coupling terms due 



to the magnetostatic fields generated by individual spin- 
waves. This coupling leads to the appearance of narrow 
instability tongues around the parametric resonance con- 
dition w(g;cos#o) ~ w o> where wo is the magnetization 
precession angular frequency. 

Spin-wave instabilities occur only for particular com- 
binations of external magnetic field and injected spin- 
polarized current. In addition, instabilities result in lim- 
ited spatial and temporal distortions which somewhat ob- 
scure but yet do not completely disrupt the precessional 
character of the original state. This robustness of excited 
precessions with respect to spin-wave instabilities has a 
precise physical origin. Indeed, the discrete nature of the 
spin-wave spectrum caused by boundary conditions in 
sub-micrometer devices reduces the number of available 
spin-wave modes which can contribute to instabilities. 
On the other hand, the strength of the magnetostatic ef- 
fects responsible for instabilities is drastically reduced, 
due to the ultra-thin nature of spin-transfer devices, and 
instability thresholds are consequently enhanced. Fi- 
nally, spin-transfer-driven precessions are characterized 
by large amplitudes and, as such, are less easily masked 
by the onset of non-uniform modes. In spin-transfer 
nano-oscillators, spin-wave instabilities are expected to 
result in increased oscillator line-widths, a conclusion 
that might explain some of the puzzling experimental 



results obtained in this area 11 1. 



To start the technical discussion, consider a ultra-thin 
disk with negligible crystal anisotropy (e.g., permalloy). 
Typically, this disk will be the so-called free layer of a 
nanopillar spin-transfer device (see inset in Fig. [TJ. The 
disk plane is parallel to the (x, y) plane and is traversed 
by a flow of electrons with spin polarization along the e z 
direction. The dimensionless equation for the dynamics 
of the normalized magnetization m(r, t) (|m| 2 = 1) in 
the disk in the presence of spin transfer is [l|, |l2( : 
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FIG. 1: (Color online) Stability diagram in (h az ,/3/a) con- 
trol plane for a ultra-thin permalloy disk. System parameters: 
a = 0.02, d = 0.6, R = 23.6, N± = 0.02 (lengths are measured 
in units of the exchange length l EX = 5.72 nm). Magnetiza- 
tion is parallel to spin-polarization in region P; anti-parallel 
to spin-polarization in region A; precessing around the spin- 
polarization axis in regions O and SW. Dashed line is an ex- 
ample of line of constant precession amplitude (cos So = 0.5) 
computed from Eq.((2|. Spin- wave instabilities occur in region 
SW. Small framed area is shown in detail in Fig. [2] Inset: 
typical geometry of a nanopillar spin-transfer device. 
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Here, the external magnetic field h az e z and the magneto- 
static field h M are measured in units of the spontaneous 
magnetization M s , time in units of (7M S ) _1 (7 is the 
absolute value of the gyromagnetic ratio), and lengths 
in units of the exchange length. The external field is 
perpendicular to the disk plane, while the spin-transfer 
torque is simply proportional to the sine of the angle be- 
tween m and e z . The parameter (3 is proportional to the 
spin-polarized current density (see [12j for the detailed 
definition) , and in typical situations it is comparable with 
the damping constant a. 

Whenever |h az - /3/a\ <N Z -N± (N z and N±_ are the 
disk demagnetizing factors, with N z + 2N± = 1), Eq.(p} 
admits time- harmonic solutions mo(i), corresponding to 
spatially uniform precession of the magnetization around 
the z-axis [l3[ (see Fig. [1} . The precession amplitude and 
angular frequency are respectively equal to: 



cosfo = 



haz - /3/a 
N x —N\ 



luq = — 
a 



(2) 



To study the stability of mo(i), consider the perturbed 
motion m(r, i) = mo(t) + 5m(r,t), with |<5m(r, t)\ <C 1. 
The corresponding magnetostatic field will be: h M (r, t) = 
—N z m 0z — N±m ± + 5h M (r, t), where 5h M represents the 



magnetostatic field generated by 5m. Since we are inter- 
ested in ultra-thin layers, we shall assume that <5m does 
not depend on z: 6m(r,t) = 5m(x,y,t). 

The perturbation 5m is orthogonal to m (t) at all 
times, since the local magnetization magnitude |m| 2 = 1 
must be preserved. Hence, it is natural to represent 5m 
in the time-dependent vector basis (ei(£), e 2 (t)) defined 
in the plane perpendicular to mo(i), with &2(t) paral- 
lel to e z x mo(t) and ei(t) such that (ei,e2, mo) form 
a right-handed orthonormal basis. The perturbation can 
be written as: 5m(r,t) — 5m\{r, t) ei(t)+5m2(r,t) e2(t). 
By linearizing Eq. ([T]) around mo (t) and averaging the lin- 
earized equation over the layer thickness, one obtains the 
following coupled differential equations in matrix form: 



1 a 
-a 1 



—N 1 - 



d_ 

dt 



5mi 
5vri2 

N ± + V\ 




1 

-1 

Smi 
5m 2 



( (s\ 
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where 



d 2 /dx 2 + d 2 /dy 2 , while (...) represents the 
z average over the thickness of the disk, and {S\ I ) 1 = 
(6\)-e 1 (t),(5\) 2 = (5\)-e2(t). 

To grasp the physical consequences of Eq. ([3j) , consider 
the plane-wave perturbation 5m(r,t) = a(i)cxp(iq • r) 
in an infinite layer (N± — 0). The corresponding magne- 
tostatic field is 1141: 



(SK) 



-s q 5m z 



(1 



1 — exp(— qd) 
qd 



(4) 

where 5m z = (5m ■ e z ) e z and 5m q = (5m ■ e q ) e g , e q be- 
ing the unit vector in the q direction. The field — s q 5m z 
is generated by the magnetic charges at the layer sur- 
face, whereas the field — (1 — s q ) 5m q is due to volume 
charges. By taking into account that Vj_<5m = — q 2 5m, 
5m z ■ ei(t) = <5misin 2 6>o, and 5m z ■ e 2 (t) = 0, one 
finds from Eq.([3]) that mrj(i) is always stable with re- 
spect to the action of exchange forces and surface mag- 
netic charges. Only volume charges can make the pre- 
cession unstable. This conclusion follows from the fact 
that the z-axis, along which the surface-charge magneto- 
static field is directed, is a symmetry axis for the prob- 
lem. Surface-charge-driven instabilities may appear in 
non-uniaxial systems. 

The two-dimensional and uniaxial character of the 
problem makes it natural to introduce polar coordinates 
(r, <f) in the disk plane, with the origin at the centre of 
the disk. The natural boundary condition in polar coor- 
dinates is d5m/dr\ r=R = 0, where R is the disk radius. 
The generic perturbation satisfying this boundary condi- 
tion consists of cylindrical spin-waves of the type: 



5m(r, cp, t) 



+00 00 

E E 

n= — 00 fe=0 



Bnk(t) Jn (q n kr) exp(inq 



(5) 
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where J n (z) is the n-th order Bessel function. The wave- 
vector amplitude q n k is identified by two subscripts be- 
cause, for each n, it must satisfy the boundary condi- 
tion dJ n (z)/dz = for z = q n kR, which has infinite 
solutions q n Q , q n i , q n 2 , . . . of increasing amplitude. The 
cylindrical spin- waves F n k(r, <fi) = J n (q n ki~) exp(in4>) are 
a complete orthogonal set of eigenfunctions of the 
operator: V j_ F nk (r, <f>) = -q nk F nk (r,<j)). The magne- 
tostatic field 5h M can be computed by applying Eq.Q 
to the plane- wave integral representation: F n k (r, <j>) — 
l/(2wi n ) Jq* exp (i q„ fe • r) exp(inip) dip, where the polar 
representation of r and (\ nk is r = (r, <fr) and q n k = 
(<lnk,' l P), respectively. By following these steps, writing 
a„ fe (i) as a nfe (i) = c„ M (i) ei(i) + c nk , 2 {t) e 2 (i), and ne- 
glecting small terms proportional to N±_ , Eq. Q is trans- 
formed into the following system of coupled equations: 



dc n 



oo A + 



+ ^ r ' P ^n+2,p(t) Cn+2,p (6) 



dt ^ A„fc 

p=0 



~T ~ ^n-2,p(*) c n-2,p ! C„i 



p=0 



C n k,2 



where: 



^4n/c — 



1 / 1 



1 + a 2 \ a 1 



-Vnk 

Vnk - K n k sin 2 9 



T^nk{t) = exp(2ioj t) - — x 



(7) 



(8) 



1 / 1 —a \ f i cos#o —1 



1 + a 2 V a 1 



— cos tfn — « cost 



^«fc-p = / rJ ri (qnkr) J n (q n ±2, P r) dr , 
Jo 



and A„ fc = r J% (q n kr) dr, v n k = q n k + i 1 ~ s nk) /2, 
K n k = — 1 + 3 (1 — s n k) /2, Snfe being the value of s g in 
Eq.gJ) for g = q nk . 

The coupling terms proportional to lZ nk (t) in Eq.© 
are the consequence of volume-charge magnetostatic ef- 
fects. They are all of the order of (1 — s q ). One has that 
(1 — s q ) -C 1 up to q ~ 1 in ultra-thin layers with d < 1 
(see Eq.Q). If one neglects these terms altogether, one 
obtains a system of fully decoupled equations for indi- 
vidual cylindrical spin-waves, characterized by the dis- 
persion relation: 



(9) 



0/(07; cos O ) = [q + 



1 - So 



(10) 



1 - 3s 



x \q + s q 



q cos 2 6 



which is obtained from Eq.([7]) in the limit a — >• 0. How- 
ever, the time-periodic coupling terms may give rise to 
parametric instabilities. Interestingly, these instabili- 
ties are governed by a small number of dominant terms, 
which can be identified by using the asymptotic formula 
Jn(z) ~ y/2/Trz cos (z — nir/2 — 7r/4) in the equation ex- 
pressing boundary conditions. One obtains the estimate 
q nk ~ it (2s + 1)/4R, where s = \n\ + 2fc. When this 
approximate expression is used for q n ±2,p in Eq.(j9|), one 
obtains: 



^ 2 : A nfc;p - ^p,fcTl i A «0;p - 1 

a = ±1 : A7u.„ = Al, = A lfe <5 pfc , (11) 



n < -2 : A± fc;p ~ A„ fc £ p , fc±1 , A ? ; 0;p ~ . 

These relations have an important physical conse- 
quence, which is best appreciated by rewriting Eq.© in 
the form: Sm = ^ m ^ s \ where: 



cW s) 



E 

|n|+2fe=s 



a n /s(i) Jn (qnkr) exp(inq 



(12) 



Under the approximation (| 1 1 [) . one finds from Eq.© 
that <5m( Sl ' is decoupled from <5m( S2 ' for any s 2 ^ Si. 
On the other hand, for each s, the (s + 1) cylindrical 
waves (namely, n — s, s — 2,...,— s + 2,— s) involved in 
<5m( s ) form a one-dimensional chain, in the sense that 
only neighboring waves in the above list are coupled. 
The absence of coupling between distinct chains would be 
complete if the approximation q nk ~ tt (2s + l)/4i? were 
exact. In that case, all the cylindrical waves in <5m( s ) 
would be characterized by exactly the same wave-vector 
amplitude. 

Instabilities are governed by the multipliers of the one- 
period map [lfl associated with the dynamics of (5m( s ). 
We have used Eqs.© and (TTT|) to make a numerical study 
of these multipliers for different chains, in order to obtain 
the instability pattern associated with each of them. The 
results for a permalloy disk with radius R = 135 nm and 
thickness d = 3.43 nm are shown in Fig. [TJ The band (O 
+ SW) in between the P and A regions is where magne- 
tization precession occurs. Spin-wave instabilities appear 
in region SW. Each chain <5m( s ) provides a distinct in- 
stability channel. In particular, chains s = 1, s = 2, 
and s = 3 (s — yields no instability at all) give rise to 
well-separated instability regions that can be neatly re- 
solved, as shown in Fig. [2ja). In general, the s-th chain 
gives rise to an instability tongue around the parametric 
resonance condition w (q; cos 6o) ~ ujq, where ui (q; cos 9o) 
is given by Eq. (fT0)) and q is of the order of the wave- 
vector amplitudes involved in the chain (see dashed lines 
in Fig. Eta))- According to parametric resonance theory, 
resonance occurs for uj — hujq/2, n — 1,2, ... . The dom- 
inant, lowest-threshold resonance occurs for n = 1, that 
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FIG. 2: (Color online) (a): Magnification of Fig. \T\ Labels 
s = 1, 2, 3 identify the perturbation chain responsible for the 
corresponding instability tongue. The pair of dashed lines 
accompanying each of the s = 2 and s — 3 tongues (one 
line only for s = 1) represents the parametric resonance con- 
dition uj (q n k', cos &o ) = wo for the largest and smallest q n fe 
in the chain. Horizontal line at f3/a = 0.15 is line along 
which the computer simulations shown in (b) and (c) were 
carried out. (b) and (c): Magnitude m^ s of average in- 
plane magnetization obtained from numerical integration of 
Eq.([T| under decreasing (b) and increasing (c) external mag- 
netic field. Dashed line represents the prediction of Eq.@ 
for sin So. Snapshots illustrate magnetization patterns ap- 
pearing just after the instability jumps. Vertical dotted lines 
are guides for the eye to compare thresholds with theoretical 
predictions obtained from (a). 

is, at to = wq/2. However, one can see from Eq.© that 
the parametric frequency is 2wo rather than uq, which 
explains why the resonance condition is uj ~ loq . This pe- 
culiarity is the consequence of the rotational invariancc 
of the problem, and is expected to disappear in situations 
with broken rotational symmetry. Interestingly, Fig. [T] 
reveals that, for a given precession amplitude cos#o (see 
dashed line), applying larger fields and currents has a 
stabilizing effect on the precession. Also, larger fields 
stabilize precessions of given frequency ujq = (3 /ex. 

To test the predictions of the theory, we have car- 



ried out computer simulations based on the numerical 
integration of Eq.([T]) by the methods discussed in Ref. 
[lij ]. Simulations were carried out by slowly varying 
the external magnetic field under constant current. As 
shown in Fig. Hfb), at large fields the magnitude m^ 8 
of the average in-plane magnetization is in full agree- 
ment with the prediction of Eq.Q for spatially uniform 
precession. Then, under decreasing field m^ s exhibits 
well-pronounced jumps, whose positions agree with the 
theoretical instability thresholds for the s = 2 and s = 3 
chains within ten percent. Beyond these jumps, non- 
uniform modes appear in the dynamics (see Fig. [2f b)), 
characterized by two-fold and three-fold patterns that are 
consistent with the symmetry of the cylindrical waves in- 
volved in the s = 2 and s = 3 chains, respectively. Agree- 
ment with the theory is also confirmed by the hysteresis 
in the instability thresholds occurring under decreasing 
or increasing external field (Fig. HJc)). 

The stability of spin-transfer-driven magnetization 
precessions has been studied in this Letter under the 
simplest conditions, namely, uniaxial symmetry and pure 
sin $o angular dependence of the spin-torque. Several fea- 
tures of physical interest have emerged: the role played 
by non-equilibrium spin-waves; the fact that instabilities 
are governed by distinct chains of magnetostatically cou- 
pled spin- waves; and the fact that excited precessions are 
not completely disrupted but only somewhat obscured 
by spin-wave instabilities. Future work will be devoted 
to extending the present approach to more general, non- 
uniaxial geometries. 
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